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Abstract 

A  procedure  for  computing  the  potential  flow  about  bodies  of 
revolution  and  symmetric  two-dimensional  forms  in  arbitrary  states  of 
motion  is  presented.  Solution  for  general  motion  is  obtained  by  super¬ 
imposing  the  solutions  of  the  boundary-value  problems  for  the  different 
components.  Technique  consists  of  formulating  each  problem  as  a 
Fredholm  integral  equation  of  the  first  kind  which  can  be  closely 
approximated,  by  means  of  a  quadrature  formula  of  moderate  order,  by 
a  set  of  linear  equations  with  a  matrix  having  a  strong  principal  di¬ 
agonal,  suitable  for  solution  either  by  elimination  or  iteration.  For 
the  general  motion  of  "bodies  of  revolution,  three  such  Integral  equa¬ 
tions  need  to  be  solved;  for  the  two-dimensional  forms,  solutions  for 
four  potential -flow  problems  are  required,  two  of  which  are  obtained 
from  integral  equations,  the  remaining  two  being  obtained  from  the  in¬ 
tegral  equation  solution  for  longitudinal  flew  by  ccnformsl  Sipping. 

Method  has  been  programmed  for  the  IBM  650,  an  automatic 
computer  of  moderate  speed  and  capacity.  The  separate  programs  for 
the  body  of  revolution  and  two-dimensional  forms  are  presented  and  the 
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application  of  computer-obtained  results  is  illustrated  vith  several 
examples* ' .A  considerably  shortened  and.  somevhat  less  accurate. procedure, 
suggested  for  use  when  an  automatic  computer  is  not  available,  is  also 
illustrated  in'  detail  and  compared  vith  computer-obtained  results. 

Introduction  •  .  , 

A  procedure  for  computing  the  potential  flow  about  bodies  of 
revolution  which  are  simultaneously  translating  at  some  angle  to  the  axlB 
of  syraaetry  and  rotating  about  a  transverse  axis  will  be  presented*  For 
a  given  body  the  method  consist's  essentially  of  solving  for  three  separ-  .* 
ate  velocity  potentials,  linear  combinations  of  which  yield  the  velocity 
potential  for  an  arbitrary  state  of  motion. 

The  numerical  procedure  for  each  of  the  aforementioned  poten¬ 
tial  flow  problems  'consists  of  computing  a  matrix  of  coefficients  for  a 
set  of  linear  equations,  the-  solution  of  which  gives  the  values  of  the 
potential  or  velocity  at  a  finite  number  of  points  along  the  body.  The' 
labor  of  setting  up  and  solving  three  sets  of  simultaneous  equations  is 
.reduced  by  the  fact  that  the  matrices  of  the  coefficients  of  the  linear 
equations  are  identical  for  rotation,  and  translation  in  a  direction 
perpendicular  to  the  axis  of  symmetry  and  that  the  operations  performed 
in  obtaining  the  matrix  for  translatory  motion  along  the  axis  of  symmetry 
can  be  utilized  in  obtaining  the  other  sets  of  linear  equations. 

A  similar  procedure  can  he  formulated  for  symmetric  two- 
dimensional  forms,  simultaneously- translating  at  an  angle  of  attack  and 
rotating  about  an  axis  perpendicular  to  the  section.  In  this  case  the 


general  problem  can  no  analyzed  Into  fovtr  fundamental  oneu,  viz,  flcrvo 
correnpondiuc  to  body  motions  parallel  and  perpendicular  to  the  plane 
of  symmetry,  rotation,  and  circulation.  By  taking  advantage  of  the 
theory  of  conformal  mapping,  the  solutions  for  two  of  these  problems 
can  be  expressed  in  terms  of  a  third,  so  tl»t  for  only  two  of  the  cases 
is  it  necessary  to  execute  the  aforementioned  numerical  operations. 

f 

It  is  well  known  that  the  Neumann  problem  of  potential  theory, 
as  the  present  potential-flow  problems  nay  be  classified,  can  be  formu¬ 
lated  in  terms  of  Fredholm  integral  equations  of  the  second  kind 
[1,  2,  3,  4,  5]."-  In  general  these  integral  equations  are  slr«ular,  eo 
that  their  application  to  obtain  numerical  solution#  baa  resulted  in 
complex  and  tedious  procedures.  Less  familiar  is  the  fact  that  these 
problems  can  also  be  expressed  in  terms  of  Fredholm  integral  equations 
of  the  first  kind  [Gj,  with  much  simpler  kernels.  Possibly  the  fact 
that  integral  equations  of  the  second  kind  possess  unique  solutions, 
while  those  of  the  first  kind  have  solutions  only  if  certain  additional 
conditions  are  satisfied,  has  been  a  strong  deterrent  to  the  applica¬ 
tion  of  the  latter  type. 

The  practicability  of  using  Integral  equations  of  the  first 
klsd  has  boon  discussed  by  the  author  in  a  previous  paper  [7).  It  was 
shown  that  even  in  oases  where  the  integral  equation -does  not  possess 
a  solution,  it  is  possible  to  satisfy  the  aquation  approximately.  A 

lfuribert  in  bracket*  designate  reference*  at  end  of  paper* 


uncLCKU>JL^L~3U 


good  example  of  this  is  von  Kerman's  method  [Qj  of  obtaining  the  ssEatly 
symmetric  flow  about  a  body  of  revolution'  by  solyinc  ait  integral  eqtto-- 
tion  of  the  first  kind  in  which  theunknu  m  function  is  the  strength  of 
an  axial  distribution  of  sources.  A  n<o+*  od  of  solving  integral  eq.ua- 
•  tions  of  the  first  kind  by  iteration  [ "v  7 j  has  the  additional  advantage 
that  the  error  diminishes  in  a  least-square  sense,  although  the  disagree¬ 
able  possibility  exists  that  beyond  some  nth  iteration  the  error  may 
accumulate  and  increase  at* some  point,  even  as  the  integral  of  the  square 
of  the  error  continues  to  diminish.  Nevertheless  one  can  obtain  a  use¬ 
ful  approximation  even  in  this  case,  since  the  errors  are  observed  at  . 

■  *  ■  *  «  «  . 

each  step  and  the  calculations  can  cease  when  the  error  exceeds  an  ac-  * 
'ceptable  value  at  any  point.  '  „  •  ■  • 

Wien  an  integral  equation  of  the  first  kind  has  a  solution, 
the  theory  [7]  indicates  that  it.  is  possible  to  find  functions  which 
satisfy  the- equation  uniformly  as  closely  as  .one  desire3.  These  func-  . 
tions  will  converge  to  the  desired  function,  however,  only  if  certain- ' 
additional  conditions -are  satisfied.  Unfortunately  the  nature  of  these 
additional  conditions  is  such  that  it  is  difficult  to  verify  in  a  prac¬ 
tical  case  whether  or  not  they  are  satisfied.  ...  ■ 

On  -the  favorable  side  is  the  fact  that  these  integral  equa¬ 
tions  are  considerably  simpler  .than  the  corresponding  ones  of  the  second 
kind  for  the  problems  with  which  we  'shall  be  concerned.  If  these  prob¬ 
lems  are  programmed  for  solution  by  a  digital  computer,  the  former  can 

be  processed  with  a  computer  of  moderate  speed  and  capacity  but  the 

&  * 

latter  requires  a  computer  of  much  greater  capacity.  With  some  sacri- 
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ficc  in  accuracy,  it  io  practicable  to  oolve  these  potential -flow  prob¬ 
lems  with  the  aid  of  a  desk  calculator  if  the  problem  io  formulated  as 
an  integral  equation  of  the  first  kind,  but  not  if  of  tho  oacond  kind. 

Experience  with  a  large  number  of  bodies  of  revolution  indi¬ 
cates  that,  for  well-rounded  bodies,  sufficiently  accurate  solutiono 
can  be  obtained  without  difficulty  by  means  of  integral  equations  of 
the  first  kind.  Tot  bodies  with  sudden  changes  in  slope  and  curvature, 
or  with  local  bumps,  the  method  of  integral  equations  of  the  second 
kind  has  been  remarkably  successful  [5]. 

The  method  employed  in  the  present  paper  io  that  of  the  inte¬ 
gral  equation  of  the  first  kind.  Five  such  equation's  will  be  derived, 
one  for  each  of  the  potential  flow  problems  mentioned  earlier.  Methods 
of  solving  those  equations  will  bo  described,  and  illustrated  by  several 
examples. 
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matrix  for  calculating  derivative  of  a  function 
semi-minor  axis  of  ellipse  or  ellipsoid,  b^  ■  b(x^) 
eccentricity  of  ellipse  or  ellipsoid,  e^  •»  o{x^) 

functions,  ^  -  f^),  gJL  «  g(x1),  hj,  -  h(x1) 

ouccessivs  anproximrrtions  to  c(x);  »  S^(x^) 
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corrected  matrices 

kernel  of  integral  equation,  1;^  ■  kfx^,  Xj) 
auxiliary  kernel,  k^  «  k'fx^,  x^) 
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auxiliary  kernel  for  ellipoc  or  ellipsoid 

unit  vector  alone  outvord  normal  to  ourfaco 

pressure  on  body  relative  to  presoure  at  infinity 

radial  distance  fro*  axis  of  body  of  revolution  to 
surface,  ■  r(xt) 

radial  distance  from  axis  for  auxiliary  ellipsoid 

radius  of  circle  into  vhich  profile  is  napped 

position  vector  of  point  on  surface  of  body  of 
revolution  relative  to  origin 

arc  length  along  body,  Manured  clockwise  Tram 
x  ■  -1  for  body  of  revolution,  counterclockwise 
fro*  x  ■  +1  for  two-dimensional  for* 

velocity  caepooentu  at  a  point  of  the  fluid 

factor  of  tangent iol  velocity  component  along 
body  for  axial  flow 

factors  of  tangential  velocity  oo*pcocnt  along 
body  far  tmhsverae  flow 

velocity  vector  at  a  point  of  ths  fluid 

mnfffiilar  Cartesian  coordinates 

seroe  of  Legendre  polynomial  of  order  M 

ordinates  of  two- dimensional  profile,  y^  «  y^x^) 

ordinates  of  auxiliary  ellipse,  y^  -  Y^i\t  *j) 

complex  variable  i  *  x  4  iy 

area  of  two  fllmeasionel  profile 

kernel  integrals  for  ellipsoid 

kernel  integrals  for  ellipse 

matrix  for  calculating  integrals  of  a  Amotion 

error  function  corresponding  to  g^x) 
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>  kernel  integrals,  1^  *  i(x^),  etc. 
h’  Jni  1  , 

scni-pcriractor  of  profile  or  meridian  section 
Legendre  polynomial  of  order  It 
first,  secondhand  rth  derivatives  of  t^( x) 
J‘G(x)ki(|,  x)dxj  Q±  »  q(xi) 
distance  from  ($ ,  0)  to  point  on  body 
weighting  factors  for  Gauss  quadrature  formula 

ffty  etc* 


*V::J  -  xi>V  etc- 


velocity  vector  of  origin  of  coordinates  attached 
to  body 

components  of  0  for  axi symmetric  body 

components  of  angular  velocity  of  body  for 
uxi symmetric  body 

resultant  velocity  on  body 

tangential  velocity  component  along  body  of 
revolution  in  axial  flow  or  two-dimensional 
form  in  longitudinal  flow 

tangential  velocity  component  foi1  transverse 
flow  along  two-dimensional  form  or  in  meridian 
plane  along  a  body  of  revolution 

tangential  velocity  component  for  transverse 
flow  normal  to  meridian  plana  along  a  body  of 
revolution 


velocity  components  of  origin  attached  to  two- 
dimensional  form 


tangential  velocity  component  along  two-dimensional 
form  due  to  circulation 
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velocities  tangent  to  circle  Cue,  rcnpcctivcly,  to  a 
velocity  of  -1  in  x-direction,  -1  in  y-dircction,  and 
unit  circulation. 

velocity  vector  of  a  point  on  (surface  of  boay 

component  of  V  along  outward  normal  to  body 

angle  of  yaw  or  attack 

arc tan  dr/dx  or  arctan  dy/dx 

Kronecker  delta,  ■  0  if  i  p  j,  «  1  if  i  «  J 

complex  variable,  i  m  f  +  iej 

imaginary  part  of  t 

angle  of  cylindrical  coordinate*  (x,  r,  o)  for  body 
of  revolution;  angle  fra*  x-axl*  to  point  on  circle 
in  conformal  mapping 

component*  of  velocity  in  circle  plane 

alternative  designation  for  x  in  Integral  equation; 
also  real  part  of 

mo*a  density  of  fluid 

velocity  potential  corresponding  to  unit  value  of 
on  auxiliary  potential  function 
factors  of  and 

factors  of  «u»d 

velocity  potential  corresponding  to  unit  circulation 
stream  function  conjugate  to  f' 
angular  velocity  of  two-  dimensional  fom 
angular  velocity  vector  of  axlsywMtric  body 
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F  circulation 

$  total  velocity  potential 

Formulation  of  Problem  for  Bodies  of  Revolution 

It  vill  be  supposed  that  the  x-axis  coincides  with  that  of  the 
body  of  revolution,  that  the  body  extends  along  its  axis  from  x  «  -1  to 
x  -  +1,  and  that  the  x,  y,  z  axes,  which  are  fixed  in  the  body,  form  a 
rectangular  right-handed  coordinate  system,  as  shown  in  Figure  1.  The 
equation  of  the  body  surface  is  expressible  in  the  form 

*  Alx)  ymA(X)AM/6  (l) 

where  r(x)  is  a  prescribed  function,  and  hence  the  outwardly  directed 
normal  vector  n  of  the  surface  has  the  components 
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/*n  -  u+r, 


where 
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Let  U  denote  the  velocity  vector  of  the  origin,  with  components 


U^,  Ug,  and  U^,  and  m  the  angular  velocity  vector  of  the  body,  with 
components  U^,  Ug,  Ug.  Then  the  velocity  of  a  joist  on  the  surface  of 
the  body  is 


7,  *  U  *3*# , 


where  v%  is  the  position  vector  of  a  point  on  the  surface,  and  the  compon¬ 
ent  of  the  velocity  along  the  normal  is 


(3) 


»  _  U,  Mrt/y  Uz  CtH.6  t04.Y  -t-  Uj  Avn,6  C6*ty 

—  ■  (  tfj  Aam,  6  -  Ufc  Uri.  Q)  (  X  Ce«t  f  +  ft.  /i£n  y)  ■ 

•  Kext  let  v{x,  y,  z)  denote  the  velocity  of  the  fluid,  vith 

components  u,  y,v,  and  §  the  velocity  potential.  Then  <j>  is  expres¬ 
sible  in  the  form'  '  * 


$  =  zufo  -  va*i  = 


(4) 


whore  <p.  is  the  potential  when  =  1,  and  =  0  for  all  values  of 
the  index  except  i  =  i,  and  the  normal  component  of  velocity  on  the 
surface  is  .  • 


it  -  2  Di  .  - 

dni  i-1  **■ 
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Hence  the  "boundary  condition  for  an  imperrcable  "boundary,  ‘ 


=.  Vs*  n 
bn  5 


yields,  by  comparing  (3)  and  (5)  and  equating  .the  coefficients  of  the 
corresponding  velocity  components, 

=  -^r,  •  =  c#i.0a>cr,  .0  tsty} 


in 
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in 
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Since  the  body  is  moving  through  otherwise  undisturbed  fluid,  we  also 


have-  the  conditions 


,$.(oo)‘=  O 


(7) 


Derivation  of  Integral  Equations  for  Bodies  of  Revolution 

Let  4-  be  another  function  harmonic  in  the  region  exterior 
to  the  given  "body  and  vanishing  at  infinity,  Ihen  Green’s  reciprocal 
"theorem  [1]  gives  ,  •  •. 

'  =  ■//•*/  & -is 


(8) 


•where  the  double  integrals  are  taken  over  the  surface  of  the  tody.  It 
vill  he  supposed  that  the  <j>!  are  selected- of  the  forma  ’ 

<t>'  «  n)t  (  i,A)  ttrt.6,  <f>'  =  <p} ix, 

=  <f>'(x,/i)  uh,9- 
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(9) 


•  A.  .  Axial  itotion 
./ 


?or  4t'  we  may  assiune  the  doublet  potential 
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which  has  the  corresponding  Stokes  stream  function 


A. 


(11) 


In  terms  of  the  cylindrical  coordinates  (x,  r,  0),  Equation  (8)  now 
becomes  . 

H' 


f*(*i !*<&) 


'  <U  =  o 


(12) 


wiidre  3  is  arc'  length  along  a  meridian  section  measured  clockwise  from 
"the  point  x  =  -1  and  2P  is  the  perimeter  of  the  section  •  But,  from  the 
relations  Between  an  axisymmetric  potential  and  its  Stokes  stream 
function,  we  have 


ji 


3?t  d-4+ 

and  hence,  integrating  hy  parts,  we'  oBtain 


At  »l‘a.  .  I %,  iST  u  =  -  f  r  g  u 
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Also  let  Ulg  Be  tlve  velocity  along  the  Body  when  the  flow-  is  made  steady 


By  superposing  a  uniform  stream  of  unit  velocity- in  the. negative  x- 
direction.  Then  •  ‘ 


34- 


34. 


(13) 


■Hence,  substituting  for  "6^/iA.  ‘in  the -previous  integral,  we  oBtain 
from  (12)  ‘  * 

P.  ....  it  ‘ 


■  .  •  {  &+**£]* 
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or  .  -  .  '  .... 

=  JP.(“V;#Wst  +  A^'d/t)  ■  (14) 

"  •  •  6  .  • 

But  since  ,  "  *  * 

‘  dtp* 

9x  tSl 

the  differential  in  the  last  integral  is  exact  and,  as  is  easily  veri- 
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•ft.1  .  A  <*-  k)  ,  ,  /  3f-  fc. 

.  -«r ^  + -gi-^. ^  =-J(r?i;  . 

Hence,-  putting  ds  =  clx  sec  jr ,  ve  obtain 

j‘  u,A***r-~  d*  =  -'2  ,'  .(15) 

-/  n* 

a  Fredholm  integral  equation  of  the  first  kind  for  the  unknown  function 

V*>;  •>■  ; 

B.  Transverse  Motion  .  ... 

The  velocity  potential  ^  satisfies-  Laplace's  equation 

a»4  ,  ,  -d*  _l_  _  >  ‘  ’  ■  • 

'3x*  +  A  <5/t  ■  .  d9a  (16) 


By  the  method  of  separation  of  variables  it  i3  found  that  <£,  may  be 


taken  of  the  form 


n>Q)  -  trt*,*)  m  6  *  (3 

where  <^.  satisfies  the  equation 

3Z(^^r)  2>Z(A  <f>r)  _  Hntr)  ■. 

■  ■  i  X 1  <>/l*  A  d/I  .  - 

which  shews  that'  r  <£r  is  a'  Stokes  stream  function. 

.  Now'  suppose  that  a  uniform  stream  of  .unit  velocity  in  the 

negative  y-direction  is  superimposed  on  the  flow.  The  flow  becomes 
steady,  and  the  velocity  at  the  body  is  tangent  to  it,  with  the 
velocity  components  in  the  directions  of  increasing  a  and  9 


jCfw"  •*  .  f  **  ,  »  M  «  j" .  f  *y  py  ^  "  T.  *•  %  m  -  *  »  *  » i*  *«  •  »  * 
,»•  >  *  ■  ‘  _  *  -  V  »  V  *«  •  •  .  %  \  . 


-1*31- 


U74  =  *24, 


U 


20 


X« 


L&lQ-  k 

<  ■*■ 

*2^  = 

*<Pr 

44 

JCiivB  ^ 

«2*  =r 

_A 

/t 

—  /iMvir 
■*  / 


) 


{!%} 


For  f  ve  may  assume  the -'doublet  potential 


A  >  _  JL  '-  AL  M4  0- 

?a  ~  r3  ** 


(03) 


in  accordance  with  (9).  Ihen,  applying  the  boundary  condition  (6), 
2<juation-(8)  yields  •.  ..*/■. 

±.(*)JA=  frJ*l 

J  A^r>nVR3;<A  J0  fi3  H  '  " 

,  substituting  for  4y  from  (18),  we  obtain 

■  ■  A;i  jf.fi] ^  . 


or 


But 


A 

drtr 


/i»  _  3  /  Al 


3  /  A)  3a 

~Uv 


(Rj)S  «V  d.Tt 

_  34  lx-  fc)  3a.  /•_£_  _  yVt  if 

ff3"  34  **/  ** 


(20) 


Tnen 


iFu  Al±  riL  \  d 4  =,  fpf  3*3(*-*>  At  _  3£*  ii  ]«/4 
J,  M.*<9  ^  an  (  *3/  L  34  3-3.  J 

'  3  J P_  jL.  J.  Zjx^)  +  Jj. 

i  **1  n  f?3  J 


_!  _  r  *(*-*)  +  A_(X~1)]P=  -  ¥ 
~  l  ft  k3  V 


«■  r. 


j 


-is-* . 


or,  more  explicitly,  from  Fq.  (20),  . 

=  -*  (^> 

*  Fredholm  integral  equation  of  the  first  kind  for  the  unknown  function 

a29  *  • ;  '  .  ■  •  • 

’’  Once  (21)-  is  solved  for  u  ,.  from  (17)  and  (18)  the  velocity 

*  bw  ,  ', 

potential  and  the  various  velocity-  components  are 


-.•  <j>2  =  4>r  <f>t  =.A(l-  *i0)  ■ 

.  .  =  [M*%?  ~  ±(a u.u)  cjty]  a** 

A  S^=  “  ( 

vu-  '  37(Aa^  ***  ... 

in  terms  of  uQ  and  duQ/dx.  The  latter  can  he  obtained  from  uQ  hy 
numerical  or  graphical  differentiation. 

.  The  corresponding  expressions  for  motion  of  the  body  in  the  . 
z-direction  can  be  derived  from  (.22)  by  replacing  6  by  0  -  if/2;  viz 

<f>f  a  <frT  MA\  6  t  St  U  -  ilto)  > 

^2  c  r  A*^r-  4=  (.«**)  my] 

J.  lha  ,(|-  U.z$)  CK&  \  (23) 

A  «  f 


U34.=  <^Y^9 


039  .  -  atfl  cwfl 
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■ C.  Rotational  Motion  ‘  ;  ■ 

The  velocity  potential  <p^  also  satisfies  (16)  and  fnay  he 
assumed  of  the  foriA  •  • 


(24). 


^  *  A<j>R  (X,  A.)  AVtuQ 
,  Let  no  assume  for  the  doublet  potential 

.  4>1  =  X  =  4 -  (25) 

*  R3  ■  R3  \  ,  l25j 

in  accordance  with  (9).-  Then,  -"applying  (6),  Equation  (8). yields 

f*  en-y  +  if^r)^ 

<■  -  *  f 

or,  from  (20), 

j\gl«‘-s<+  <»> 

A  •  _| 


a  Fredholm  integral  equation  of  the  first  hind  for  $  .  For  elongated 
bodies  a  preferable  form  of  the  equation  is  obtained  from  (26)  by.  ‘ * 
■writing 


_  r  3X>l!  .  r  3X(X-t)rt.*  /t*  <1  1  J, 

Ji  f~  Rs  t  •  aft*  J 


rjf  ^ ~  ■£  Hy ~  \  £r *y *.  "  r  *  ■*  *'*  '^T  *  4  *~r *£ ' 

1*1  -  '  ’ v  T/.  ."7-  /■'  .  *  "**  *  •  *'■  ‘  ’  y  *  *  . 

1 ' .  -X  w\ ;v '. ■■>  ti£Em> -^-y-y>  w.  -v* yy-l 


But  from  the-  identity 


J[  $  - *l&isi&>£. ]JW 


vc  obtain. 


-,.t  R/  1  .iflM  jr  ]  ■» 


Hence,,  from  (27)  and  (28),  we  "obtain 


jW.-f*.)  ^  [  R*-3kS  %«-t)  18M  t-  f  '  ^  ^  (29) 

•,.|  *  f  fif  1*  ^  4X.  J.  *  -l  RJ  h 

the  desired  integral  equation,  in- which  *<$>+*  is  the  unknown  function. 

The  velocity  potential  and  the  velocity  components  are  then  given  by 


the  expressions  ' " 

<f>^  =t  A.  4*  A* 6  •» 

^  •  j  .  ( 

jl  fit-.  =  4:  ■  .  J 

A  *e  •  *  .. 

The  corresponding  expressions  for  rotation  of  the  body  about  the  z- 
axi3  cad  be  derived  from  (30)  by  substituting  0.  for  8;  viz 


<f>t  =  -  Ji<PA  tCrC 9 

TZ*  - 

±  iii  =  A  ^  0 

A  3fl  *« 


Formulation  of  Problem  for  Symmetric  Two-Dimensional  Forms 


Only  the  case  of  a  two-dimensional  form  symmetric  with  respect 


to  its  chord  will  ho  considered.  The  x-axis  will  'be  token  to  coincide 
vith  the  chord  and  the  y-axls  perpendicular  to  it  at  its  midpoint', 
with  the  chord  extending,  from  x  «-l  to  x'=?  +1.  The  ecjuatlon  of;  the 
upper  half  of  the  profile  will  he  expressed  in  the  form  y  ■=  y(x), 
where' y(x)  is.  a  prescribed  function.  -Then  the  outwardly  directed  ■ 
normal  to  the  contour  has  the  direction  cosines  (-  sin  y,  cosy) 
where  y  =■  arctan  dy/dx.  Also  let  s  denote  arc  length  along  the  upper 
half  of  the  profile,  measured  counterclockwise  from  the  point  x  *  +1, 
and  3P  the  perimeter  of  the  profile.  Then 

•  '  .52  =  —  =  —  Q-tmy 

■  >  ZA.  bit  '  ' 


-  —  =.  — 
iA,  »«■  * 


See  Figure 


*  Let  U,  V  denote  ■  the  ■  veloci  by  components  of  the  origin  of  the  ■ 
body  and  to  the  angular  velocity  about  an  axis  perpendicular  to  the 
aection.  It  will  also  be  assumed  that  there  is  a  circulation  V  about 
the  body.  .Then  the  components  of  the  velocity  of  a  point  on  the  contour 
are  U  -  u3y  and  V  +  cox  ,  and  the  component  o‘f  the  velocity  along  the 


normal  is 


=  -  UMrtiy  +  +•  u)  (x  ttriilf  +  ^  (32) 


Note  that  this-  convention  is  opposite  to  that  used  previously  for  the 
body  of  revolution. 


>y>p  runr-n  *'^j 


Mi, 

Sn 
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The  velocity  potential  |>  is  expressible  in  the  form 


§  =  u$,  +  V4*  +  r<*. 


‘  (33) 


Then  the.  boundary  condition 


M  &  v. 


£>.TV 


(di) 


yields  the  separate  boundary  conditions  . 

Mi  _  <M±  _ 


=  >  p^  =  C**F/'  =  %  frxt-Y-h  yAlrtiy . 

*  3  ji  ^  *i 


(3S) 


in. 


Derivation  of  Integral  Equations  for  Two-Dimensional  Forms 

It  is  proposed  to  solve  these  potential  flow  problems  by 
determining  <f>,  and  <p.  from  Fredholm  integral  equations  of  the  first 
kind,  and  then  obtaining  <f>2  ,  and  ^  from.  <pt  .by  means  of  some  conformal 
mapping. theory.  . 

A.  Longitudinal  Motion  ■  •' 

As  before,  we  begin  with  Green’s  reciprocal  theorem,  but  in 
its  two-dimensional  form,  . 


=jV  | u. 


L* 


in 


(3B) 


where  <j>'  is  another  function  harmonic  in  the  region  exterior  to  the 
contour  and  vanishing  at  infinity.  When  <f>  a  <£>;  wo  may  assume  for  <p( 
and  its  corresponding  stream  function  ^/rthe  doublet  Junctions  . 


i 


H 


p 

!v 


m 


W 
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w 

f  f 


>V 

H  * 

y, 

t 


-fe 


rm 


3TO\.uX.VUtma, 43B&3&^^TmXXSZ-XjWW^ 
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4,'=-^,  r-r-JL'y  »*=(«- «Vlf*. 


(37) 


h*  '  T  /** 

Then,  substituting  for  dfy/Jri  from  (35),  applyiTig  the  Cauchy-Riemann 
equation  '  .  - 


in 


"irf/' 


and  integrating  the  left  member  of  (36)  by  parts,  we  obtain 

-f*  2L  d5  js  —  ip'jlun  y  ds  z=  f  (f>/ ds 

J.  3*  H*  J.  r  J(t  3<J 


is  R* 


How  let  U  denote  the  velocity  along  the  body  when  the  flow  is  made 
steady  by  superimposing  a  uniform  stream  of  unit  velocity  in  the 
negative  x-diirection.  Then  ' 


]L£l  -  U1(f  +• 
is  •  ** 


•and  hence  ve  obtain  from  the  previous  equation 

z-i 


p  p 

—  J.d  ( (Vutdnv 


or 


-/ 


(38) 


a  Fredholm  integral  equation  of  the  first  kind  for  U., . 

»  J5 


E**? -jr?-  h  m  *  f  /  '  »■  r*r  -  „*  ■*  * 

?>'/>■& ■  Y* 


T-v 


*-.1  a-S.  —  V.JTk-A  m*A**J\.  1,1  ^.'l^Ut.J 


jg=f 

-  if-W* 
f  kesSJS 

-  i 

t  1  ._ 

,«arr*' .  -^sb^Vj^www 
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B.  Rotational  Motion 

Ve  take  ♦  =>  and  assume  for  $'  in  (36)  the  doublet  function 
and  its  corresponding  stream  function  V  • 


'  3  ■»;+’-  ^ 


X 


*2 


(39) 


Then,  applying  the  boundary  condition  (35)  and  the  relation 


l±'fc  ££' 

fin  *■& 


Me  obtain 


<54. 


Ja. 


- jf> ♦ * *8*]*  «»> 


a  Fredholm  integral  equation  of  the.  first  kind  for  ( )  • 

For  elongated  forms  a  more  suitable- equation  is  obtained  from 
(48)  by  writing 


(41) 


*i  -  s  ■ ►  ^ *•  • -z 


and  consequently 


V 


’I 


&t2tt^s2z*i’&xxu  i 


-,a^uv.«  .  'tt,??**. 


Applying  this  result,  (4l), becomes 


Sty-  =  *$  {42 
the  desired  integral  equation,  with  (.-£.-  I ).  as  the  unknown  function 
Solutions  for  Transverse  Motion  and  Circulation  -  Kutta  Condition 


The  given  profile  can  be  mapped  into  a  circle  of  radius  ro 
about  the  origin  by  means  of  the  transformation  . 

-  •  -(43) 

where  z  =  x  +  iy  and  £  =  £+.ii£  are  complex  valuables  in  the  planes 
qf  the  profile  and  the  circle,  respectively.  Then,  denoting  the  com¬ 
ponents  of  the  velocity  in  the  £-plane  by  (/< ,  v),  the  relation 
between,  the  complex  velocities  in  the  z-  and  ,$ -planes  is 

•  ii-.iV=  (fi-iv)  ‘  (44) 


Let  us  again  consider  the  case  of  the  body  at  rest  in  a 
uniform  stream  of  unit  velocity  in  the  negative  x-direction.  Then 
•since  z.  =  £  and  d^/d£  =  1  at  infinity,  it  is  seen  from  ( 44)  that 
the  flow  in  the  }  -plane  about  the  circle  i3  also  due  to  a  uniform 
stream  of  unit  velocity  in  the  negative  ^-direction.  Henco  the 
tangential  velocity  on  the  circle  is 


U)c  -  2.  amv  Q 


(45) 


where  ©  is  the  polar  angle  measured  counterclockwise  from  the  J-axis 


>_v-a  v^DaaSsSST.aaas: 


^^C&aKtik'a&iz  t£«&*  '  ivi-V 


Hence,  taking  the  absolute  values  of  the  members  of  (id),  we  obtain 


U|s  =s  2 H0Avn,Q. 


which  gives 


H-K0  -  i  U,sds  ~  f  U,suU£.y.dx 


2 A0(i+tUrt6)x  J.Ulsd s  =  J  UIS  ^r-4x 


'.  ,  ,  ztie  (  i  -  otty  =;  /*  u)sJs  =  J  OlsMcr-dx 


■  Thus-  (4?>  gives'  tq  and  9  as  functions  of  x. 


How  consider  the  case  when  the  body  is  at  rest  in  a  uniform 


stream  of  unit  velocity  ih  the.  negative  y-direction.  The  tangential 
velocity  on  the  circle  is  Ugc  =*  -2  cos  0,  and  (44)  and  (46)  now 
give  for  the  tangential  velocity  on  the  profile 


U2S  =  -  2Aa  cone  •  4^  -  -  U,^C4rt  6 


"  Thus  the  solution  for  transverse  flow  is  given  by  (47)  and  (48)  in 


terms  of  that  for'  longitudinal  flow. 


The  result  when  the  flow  about  the  body  id  .due  to  unit  '  • 


circulation  is  equally  simple.  The  circulation  about  the  circle  is 


also  unity,,  so  that  the  tangential  velocity  about  the  circle  is 


given  by  27TrQ  Uyc  =  1.  Hence,  from  (44)  and  (46),  the  tangential 


velocity  on  the  profile  is 


it  _  j_  M  =  U,s 

VS  XI T  Js  HD' /igAim-d 
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For  profiles  vith  a  sharp  trailing  edge  the  circulation  T 
con  130  determined  from  the  Kutta  condition  that  the  resultant  velocity 
at  the  trailinc  edge  due  to  transverse  flow  and  circulation  should  be- 
finite.  From ‘(48)- and  (49)  the  resultant  velocity  corresponding  to 
©.  «  TT  is  seen  to  he  ,  . 


vuZi  +  rovs 


'  '(  ZV/t°  +-  "zf) 


.But  at  a  sliarp  trailing  edge,  d9/ds  is  infinite  and  hence,  in  order 
to  satisfy  the  Kutta  condition,  we  milst,  have 


T  -  -  ‘nrn0V  .  (50)  . 

Hence,  from  (47),  (48)  and  (49),  if  such  a  profile  is  in  translational 
motion,  the  resultant  flow  velocity  UR  on  the  'surface  relative  to  - 
the  body- is  .  r 

uR  =  uuls  +  vuzs  +  ruV5  =  uls  [  u  -  v  ( ***  +  ^&)j 
=  u,s[u*v/i5grj  '  -  •  • 

=  uls  [I^'s  AXC?J*- 
JjJ  Vis  •*“/  _ 


in  which  the  i-  sign -refers  to  the,  lower  half  of  the  profile  and  the 
-  sign  to  the  upper  half.  .  ■  ... 


:V 
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Procedures  for  Solving  Integral  Equations 

•  We  arc  concerned  vith  solving  integral  equations  of  the  form 


J*  JU)JX  =  -f  (^) 


(51) 


in  which  k(£  ,  x)  and  f  ( J  )  are  given  and  g(x)  is  the  unknown  function. 
If  the  Integral  in  (51)  can  "be  expressed  as  a  finite  sum  by  means  of 
some  quadrature  formula,  it  would  assume  the  form 


f  KjX j9j  =  *i 


(52) 


where  the  R.  are  weighting  factors,  g  —  g(x  ),  k  =.l(x.,  xj,  and 
v  J  J  ij  *  i  0 

f^  =  f(x^) .  .  This. would  reduce  the  problem  to  that  of  solving  a  set  of 


linear  equations  for  the  unknown  ordinates ,g..  • 

„  J 


The  kernel  f unctions, k(  x)  in  the  integral  equations  that 

we  presently  wish  to  solve  all  have  the  property  that,  for  the  elongated 
bodies  of  most  practical  interest,  they  are  small  except  in  the  neigh¬ 
borhood  of  x  ='£ .  .  Although  this  property  is  desirable  for  obtaining-' 
approximate  solutions  of  the  integral  equations,  as  will  be  shown.  It 
appears  to  make  it  necessary  to  use  a  quadrature  formula  of  very  high 
order  in  order  to  obtain  an  accurate  representation  of  the  integral. 

An  essential  feature  of  the  proposed  procedure  is  that  it  is  possible 
to  modify  these  equations  so  as  to  avoid  this  characteristic  of  the 
kernels,  so  that  the  resulting  integrals  can  "be  represented  accurately 
by  quadrature  formulas  of  much  lower,  order.  For  this  purpose,  (5l)  is 
written  in  the  form 


"  .  sfC*) 

z(t)  - 


v, 


N 


■ '  *0.'.  .'.; — .*  ..  j  ■•'■»?C1.-. , 


‘  .^..  .^^  ■t.-i.-w,1.  ...  ... — ;»*,■  *-i'>-ai^/_'_i^jv.-.i..>-i>.  jx.  -aLs 


**.  c.-:u.  ^ _ ■_ 
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■<? 


vhcre  g'(|)  denotes  the  derivative  of  g(£).  Wd  will  assume  that,  the 
derivative  at  {  »  x^  can  he  expressed  .in  terms  of  the  valuee  of- 
the  function  g^,  g 2>  ...  ^  at  x^,  x^»  ...  by  means  cf  a  numerical 


differentiation  formula  of  the  form  '■ 


?t  =  £  a«/»j  ■■ . , 


4  m) 


Since  the  integrand  of  the  first  integral  in  (53)  vanishes  ■ 
when  x  ■  |,  its  peaking  property  has  been  eliminated,  although  It 
appears  that  the  difficulty  has  simply  been  transferred  to  the  other 
integrals  occurring  in  the  equation.  If  it  is  necessary  to  employ  a 
quadrature  formula  to  evaluate  these  latter  integrals,  this  will  be 
accomplished  by  writing 


jib*  (x- HLk! *>*) - Wwi j*  +  J'(V 


(55) 


where  k’(£,x)  is  an  auxiliary  matrix- to  be  selected  such  that 
k(x,x)  •»  k’(x,x),  and  I‘())  «  j|  k’(£,x)dx  and  j'(J)  ■  J*  (x-  J)k’(),x)dx 
are  integrable  in  terms  of  known  functions.  It  is  clear  then  that  It 
would  be  possible  to  evaluate  the  first  integrals  on  the  right  of  (55) 
by  a  quadrature  formula  of  low  order.  It  .-.rill  be  seen  that  the  inte¬ 
grals  of  two  of  the  four  kernels  derived  above  can  be  expressed  in 
terms  of  the  integrals  of  the  other  two, .  so  that  only  two  such  auxiliary 
keraelB  k’(  £  ,x)  need  to  be  found. 

Gauss's  quadrature  formula  is  a  convenient  and  accurate 


-TT-T  HXXJOXSJt  XJI  XJIX2\X£IXSZ 


>4T  at,  *f~\+Sl+*  .  4. . 


aJ  I 


method  of  evaluating  these  integrals.  The  formula  way  be  expressed 
•in  the  form  .  -  ■  • 


...  i  i 

- 1  I 

OJ  i 
‘■J  I 
;*’1  . 


li 


’  .  .  •;  ?Ft*U*  =  ■  -  *  (») 

.  ■ .  j*  j«/  .  • 

where  the  x.^  are  the.- zeros  of -the  Legendre  polynomials  of  degree  ? 
and  the  are  weighting  factors,  the  values  of  Xj  and  It^  have  been 
tabulated  [9]  for. values  of  N  ftora  1  to  16.  The  values  -for  N  =  11 


and  16  are  given  in  Table  1. 


Table  1 


Abscissae  and  weighting  factors  for  Gauss’s  quadrature  formula 


-0.97822366 
-0; 88706260 
-0i7301.5201 
-0.51909513 
-0.26954.516  ' 
0 

0.26954516 

0.51909613 

0.73015201- 

0.88706230 

0.97822836 


.  0.05566857. 
0.12558037 
0.18629021 
0.23319376 
0.26280454 
b.  27292509 
0.26280454 
0.23319376 
0.18629021 
0% 12558037  • 
0.05566857 


-0.98940093 

-0.94457502 

-0.86563120 

-0.75540441 

-0.61787624 

-0.45801678 

-0.28160355 

-0.09501251 

'0.09501251 

0.28160355 

0.45801678 

0.61787624 

0.75540441 

0.86563120 

0.94457502 

0.98940093 
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0.02715246 

0.06225352 

0.09515851. 

0.12462897 

0.14959599 

0.16915652' 

0.18260342 

0.18945061 

0.18945061 

0.18260342. 

0,16915652 

0.14959599 

0.12462897 

0.09515851 

0.06225352  j 

0.02715246  I 
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Applying  Gauss’ b  quadrature  formula'  and  (54)  and  (55)  in  (53),  we 


obtain 


ft,  h ij  ~  ^  >■  =  +  s,i  I**- SCI  +tt/j  tJrT0  ^ 

=  O  if  itj  ,  =  >  ~  j  i  X,  =  !<*/),  X,  *  Jt*i) 

Sj  *=  s  Rj  fcij ,  *n  *■  •?  rj  - 


l>.L~nJ.~-a.jCiiXL>;U>^v>^cX.L^  LXUyfa,X\:aL5i.L>. 


I*.*^*^  — *  WOiV»*»  «>•  w.>* A  *a£C *•£'.££.  Wt >W«> ^  Mv»fri  tu vA  — .'  '•  4,>>  *  J 


6 


If  the  quadrature  formula  is  of  low  order,  thi6  set  of_  equations  can 

he  solved  exactly  for  the  values  of  g..  It. appears  to  he  preferable, 

J. 

however,  for  reasons  already  discussed  in  the  Introduction,  to  solve 
integral'  equations  of  the  firBt  kind  by  successive  .approximations,  as 
will  be  illustrated  in  a  subsequent  section.  _  , 

■An  approximate  solution  of  (51)  may  .be  SerlvfeA  as  ■follows. 

’  Since  k(|,x)  peaks  sharply  in  the  neighborhood  of  x  »\j  and  is  gs»11 
elsewhere,  and  g(x)  varies  very  little  from  g(|)'in  this  region  in 
which  the  value  of  the  'integral  is  principally  determined,  only  a 
small  error  would  be  introduced  if  g(x)  were  replaced  by  g(j )  in  (5l).. 
trills  gives  ,  ' 

'  .  '  \  9<l)  £,/»(*,  *)'!*  fU).  •  /  ■_  ‘ 

whence  we  obtain  thb  first  approximation  . 

•  ,  »<*,  =  }j±  ;  '•  ‘  (59) 

1  ■  r  cx)  .  •  ,  ■  ' 

■  We  can  now  define  a  first  error  function  . 

E,<?)  =  I(T)  ifft)  "  J->  l'(x)k(t‘x-,</yJ'  •  ■  '  ' 

which,  by  (5l)  may  also  be  written  in  the  form 

.  *  »•  *  *  ■  ,  ’ 

•  -  E,t|)  -  — .  (60) 

Since  the  kernel  is  the  same  as  in  (5l),  an  approximate  solution  of  (60), 
Indicated  by  tile  same  argument  as  was  used  to  Justify  (59),  is 

G2(x)  -  B^x)  ■=  EjW 

or,  from  (60),  "  ■  .  . 

=  9,c*i  *T7r,  -Jj 


.Ei.AbX£3tEjiu^\j&u«  i^. 


»#-^Jf^*«^i**,W*(  ,J^*VV  *1 


<*  *JU.*.'^* /.UlOtfc'MPrX 


The  foregoing  process  can  'be  continued  successively  if  ve  define 

.  E,«)=j^5  I  3.  (j>  -/'  9,.<«  k  (  i.'J  <**r  .  (f 

Vhlch  may  be  written  as  the  integral  equation  .  , 

£7F)  (.i,C*tx)i 


an  approximate  solution  of  which  is 


W(,°  -  $*’<*>  “  *M" 


(62) 


Hence  ve  obtain 


•  1  M  58  t* ( ^  +  r(l)  f  f ( ^  <x?  klit  ’.  (63) 


°r>  by  (57), 


■  u>/s-u  : 


the  desired  iteration  formula.  This  last  result  may  also  be  inter¬ 
preted  as  an  iteration  formula  for  solving  the  set  of  linear  equations 
(57).  In  solving  a  problem  the  iteration  should  certainly  cease  when 
the  .  of  largest  magnitude  begins  to  increase  with  n. 

,  It  is  now  clear  from  , (57)  that  in  order  to  solve  an  integral 
equation  of  the  first  kind  by  the ' iterative  procedure  of  (64)  it  is 
first  necessary  to  evaluate  1^  and  J^.  Methods-  of  accomplishing  this- 
for  the  five  integral  equations-  derived  above  are  treated  in  the 
following  sections. 
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A.  Axiuyr.mctric  flow  about  a  body  of  revolution 


‘The  velocity  along  the  body  is  to  be  found  as  the  solution 


of  (is)  in  "which  -we  way  take 


<j  cc)  „  -  u,j  A**r ,  -*<*)  -  2>  fe'-’f i,  *)'=  A*(*W-V'+A*(pf’A( 65\ 


In  order  to  express  (15)  in  the  form  (64)  it  is  necessary  to  evaluate 


the  quantities 


>  x  {$)  =  i  Kw(t,*)<«x  ,  '  .  j'u-4)hw(t;^«l* 


This  will  be  done  in. two  different  ways. 


The  basis  of  the  first  method  is  the  observation  that  the 


exact  solution  of  (65)  for  a  prolate  spheroid  of  eccentricity  e  is 


<}(X)  =  ~  5  A  =  p  j><  *  -££]  ^  b  »  f~- t*  (66) 


This.  suggest s  as  .the  choice  of  the  auxiliary  kernel 


MM*  <U*-tf+A?]~vtr  :  K-'^ujrjz 


*  ':*  *<  t)  lllL 


Here  rB(  £,x)  represents  the  prolate  spheroid,  the  major  axis  of  which 


extends  from  x  =<  -1  to  x  *>  +1  and  which  intersects  the  given  body  at 


x  *  J>  This  spheroid  has  the  eccentricity 


e  =  i//-  fc*  ,.  '  b  = 


Wc  obtain  then  from  (65)  and  ( 66)  - 

I£L-j'kE(*j,*)Jx«  At  ,  At*A[e (*;)] 


Furthermore,  by  direct  integration,  we  can  obtain 


Jel  =  H*-*0  he(Xi,x)Jx  =  x •  B{. 


te— 1^;  s* 


t  «S*jc-v*  fc — Mj&ttVpVaWttA  ;-v->Cy.- 
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Then  1^  and  J.^  cun  be  obtained  frcm  (55)* 


The  second  method- consists  of  employing  the  transpose  of  the 


kernel  K( j,x)  as  the  auxiliary  kernel: -viz- 


,-■>/« 


Then 


’  /->■  .  1+1 


i/  *  C *•'<*« x)it  ~ /jslijj*'*  [fvijit&tm* 
j/  =  f'(*-n  k'a^ *  -  ^ l  jjjTuh^l^  ~  J  •  (?3)  .  •  • 


(71) 


(72) 


vhence  1  end  ^  can  be  obtained.  The  latter  method  has  -the  advantage 
of  avoiding  the  calculation  of  an  additional  matrix,  but  it  is  less.  , 
accurate  than  the  former  method  near  the -ends ’of  the  body.  > 

B .  Transverse  flow  about  a- body  of  revolution 

Let  us  take  g(x)  -  ^  and  f(x)  =  4  in  W  w  need' 

to  evaluate  ‘ 

i  (t)  -  j1  + 


7tH) 


-»  Hs 

J<  Cx-Da1?*)  I'M2-- }/!.*'+  ' 


vhere  G(x)  |j  (r2)-  But  from  the  identic 

_ _ _ _  ..  _  T3  jxL  n  *  J 


~2  dx 

Kraa,*)  - 


ve  obtain 


i,c«  =  i.®-* 


M&&  *?mm  ,«,-4»«. '  »-»'.  ^jifeOi’U  ,» .$.\  »  *rffc»WJU' %**.!»«,* 


*  A'  ^ r  m 
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Furthermore,  subtracting  (26)  from  (29)  yields 

X  xk^lXl,*)**  -  J  xkM(.*i,X)  dx  ~  ‘fXi  -  &j 
*  “*  ,  „  “* 
where,  frorf  (57), 

1  .<»>, 


0) 


at  =  j  h”1  c*..,*)  &<»)  a*=z  hij  g-j  . 

Since  also,  from  (-74),  .  -  • 

.  ■  '  '\  S'  *i  ( *«•  >  *)**  -  j'  x,  k0)( x,., *) dx  ~  v %: 


.(75) 


Vo  obtain 


^21  =  Jtl  ~  Q». 


(76) 


Tim's,  by  (74),  ' (75)  and  (76),  I  and  J  are  expressed  very  aiaply  in 


terras  of  1^  and  J^. 


It  is  seen  fron  (22)  that  it  is  necessary  to  evaluate  d(ru,,e)/dx 


In  order  to  determine  the  tangential  components  of  the  velocity.  Although 
this  derivative  can  be  'obtained  graphically  from  the  values  of  Ugg(x),  a 
numerical  differentiation  formula,  derived  in  Appendix  .1,  is  available  .- 
when  the  values  of  u2Q  are  given  at  the  abscissae'  of  the  Qauss  quadrature 


formula.  For  example,  the  expression  for  U2q  in  (22),  in  terms  of  the 


matrix  tabulated  for  N  *  1&  in  Appendix  1,  is 


~  Uas  <*i)  =  “  Yi  «,y  Aj  U20(Xj) 


(77) 


C.  Rotational  motion  of  a  body  of.  revolution 
In  (29)  we  may  take  , .  ’ 

=•$*+*,  Ri)=-*t'-£2GMft°}(L*)?l* 

Then  the  kernel,  and  consequently  the  corresponding  matrix  h^j  are  the 
same  as  for  transverse  flow.  It  remains  only  to  evaluate  f^  ■  f(x^) 


v.r  :_v..r  . 
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I,  ^ .  and  can  then  ho  obtained  iron  (55)-.  '*’..  .  ■ 

Xn  the  alternative  method-  the  auxiliary  kernel  ic  the  transpose 
of  the  Given  one*  Then  . 

■'  x«,~  £  '1 it  +  "*V.  if'  .<«}  • 


*  J.'  .»<»->■•)  J  U,  i- ,  (’rtiftlL  .  ’  • 


whence  1  and  can  be  obtained.  As  in  the  axisyinmetric  cace,  the 
latter  method  avoids  the  calculation  of  the  additional-  matrix,  but  it 
is  less  accurate  than  the  former  method  near  the  endo  of  the  profile. 


•  . flotation  of  a  symmetrical  two-dimferioional  form 
In  (42)  we  may  take  ■ 

.  4( £)  -.  tr\+  zjr  G (vU* 

f  %y  >  ■ 

Crti)'  s  ?  3*  >’  kU>fxJt*)  c  i 


Then  subtracting  (42)  from  (40)  giveB 
I 


XtL  =  J' k,*/(xi,x)Jx  * -  irXi  +  jJ  (x~  t)  J*' 

or,  from  (5?)*  , 

*»,  «  ii-Oi  *  *  «*>*)  *  fl.  jU»%  |M)- 


Also  we  havO  •  * 


«/ 


=  /*£**  „  ^£  4  x,u  iriu  * **jQ  j  Jx 

-»  *•  '  Rx  -  R1  *  **■  R3-  J 


'V-'.'.V'V-','-'  s  .v  ,v.-.  ..  .....  ~T*.  .  . 
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zas:  sj?  jumkn-v  a^.aw.a.«'?ij>x\  uk -a  k  ji.yx  ux.  ,  :\l^*  ^.h  _^\ 


*W  i  **■>■*>«*•»■*,  *«*r-  >  ,-MX  *•*»**>**> av -^.a*  rtv*.i«j  *»  • 


'  1  ■  |  •  • 

ft  }■**?')  J*  =  xyl^  **  ,  St  ><,j*-«  ^  =  £4, 


Hence 


JU  =  *»  (  J>1  “  5^  +  T°“  2  h‘i  [*?;+  (*i~  *0  &)  ]  .  (86) 


Lastly  it  is  seen-  that  f(£)  is  also  expressible  in  terms  of  Q  ( }). 


We  have ' 


jc  TTX,  +-Z  Qi 


’  (87) 


The  tangential  component  of  the  velocity  is  now  given  hy 


(~f).  =-^.2^  ^O+Sj) 


Summary  .of.  Procedures  for  Solving  Integral  Equations 


1.  General  Procedure 


..  klxt  ,-x;  <f»  k-  ^  h,-}  $j  *  A 

-  S*i  +  5- ( *»  -  f  bH  3"^»  =  rT 

.  h,j  »  txs~i{)  +  a,J  ^J«*“T') 

x,  *  i;  -  £■<*-*)*<*», »*d* 

-'■  Si  »  2  >  ri  “  f  *4  (ltrx0  *>«i 

.  .j^  %'x^.di^ j 


.  — ■_  4!  >.  -  ■_  <_%VV\.VLVr.^aiVu%7iA.XVu^<\flru%XVrjV\%^*«’l  VLiVL*. 


2<  BodlC3'  of  Revolution 

.  a)  Axial  flow  ,  '  "•.•••• 

k%  =  >  *4  =  f*^;)***/  ;  f{  =  2.'  '- 

X'ti  -  5a  —  Ai~  5ft  >  J/i  —  f»£  —  *£  ®i. -  T^x  -  ■ 

Al“  e7a.[2ei~  b,*-tn  .’\£j.J->"  Br^  ^£tci  "  tz+ J 

•  ■  ■  =  I  -  e:*  «  — ■ 

1  l~X,z 

Sei  =  2,Rj  KEij-  TPx  =  fKj(Xj-*i)  keij-  - 

■v  v  __  b,*(i-*^  ■  .■  ;  ' 

■  "  ei5'  {(x--xpa+b{t<.-*j?JJ/2  ■ '  • 

b)  Transverse  flow 

.  •■  .  5.  .-  4  -*.  '  i:  ■••$  •-  '  i  •' 


ai  -  f  <% 


c)  •  .Rotation  -  ^  same  as  for  transverse  flow. 

f ^  =  4x^  +  . 


3.  Symmetrical  'iVorBimensional  Forms' 


a)  Longitudinal  flow 


fc"  -  a?  •  -  tx‘'v *  >t  ■■  f'“T 


T  .  -  = 


0,-5.!.  J”.i  “  T.  ■  XiD,'-Tgx 


F9 


1 


1 
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b)  Rotation  .  •.■■;•* 

'  *<*>  *  T%--'  ,  k%  “  [  **i  ~  *»*"+  *  lfi~t0  *** 

^  =  7T fi  +  zti,  +*f  ’  &<*)  ?.  *  *&  ; 

■  **•  ,  •  •  • 
•  ..  xtl  =  j;t  -  Qi  +  *i  (•r«i“ir9  -  • .  • 

Jti  =  . 

»  i  *  ** «  ' 

•  *  •  *•  *  ,  ,  « 

*  ,  *  % 

'4.  Short- cut  procedures  .  1  • 

If  an  automatic  digital  computer  is  not  available, 
sufficient  accuracy  for  most  purposes  may  be  obtained  with,  the  ■  •  .'  . 

Xolicving  briefer  procedures: 

a)  T;&e  fox*  the  corrected  matrix  •  .  . 

'  :  iiij  =  Kik,JfT,jir.--si)  .  ■.  : 

Consequently  it  woyld  not  be  necessary  to  compute  the*  and  for  any 
of  the  integral  equations.  "  .  ’■ 

b)  Use  the  transpose  of  the  matrix  k^j  ,as  the  auxiliary  matrix 

to  obtain  1^  by  means  of  (72),  except  x’or  i  =  1  or  N.  *  Obtain  1^  ^  and 
1^  g  by  means  of  the  more  accurate  procedure,  using  the  ellipse  auxiliary 
matrix,  summarized  .above.  * 

c)  T.ike  N  =  11,  rather  than ’the  value  H  =  16  which  is  being  used 
in  programming  the  problem  for  an  automatic  computer. 

It  is  estimated  that,  with  these  short  cuts,  the  resulting 

■  labor  is  only  one  fifth  of  that  of  the  complete  procedure.  The  consequent 

*  »A 

loss  in  accuracy  can  be  judged  from  an  illustrative  problem  which  Is 
solved  in  both  ways  in  the  following  section. 
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Appli cations  '  ! 

Ere  acthod  of  sblving  the  various  integral  equations  encountered 
above,  applying  the  Gauss  16-point  quadrature  formula  with  the  ’ellipse^ 
correction,  and  solving  the  resulting  sets  of  linear  equations  by  iteration, 
has  been  programed  for  the  IBM  650  computer.  Two  programs,  one  for  bodies 
of  revolution,  the  other  for  symmetric  profiles,  are  presented  in  Appendix  2. 


A  test  of  the  procedures  was  given  by -application  to  the  form 


=  0.04(l-x4j  or  y^  a  Q.04(l-x4) 


■which  is  shown  in  Fig.  3.  Because  of  the  limited  capacity  of  -the  computer, 
certain  preliminary  calculations  are  made  by  hand,  as  is  illustrated  in 
Table  2  for  this  form.  These  data  are  then  introduced  into  the  computer 
which  operates  on  them  to'  obtain  the  various  mtrices  for  the  different 
motions  and  solyes  the  resulting  sets  of  linear  equations.  The  solutions  • 
arc  given  in  Table  .3.  ...  ‘  ' 

The  subsequent  calculations  for  obtaining  the  solutions  for  :  * 
transverse  motion  add  circulation  for  a  two-dimensional  form  and  for 
combining  solutions  and  obtaining  pressure  distributions  have  not  been  • 
programmed  for  this  computer.  For  example,  the  pressure  distribution  on 
a  body  of  revolution  in  steady  motion  at  an  angle  of  yaw- rf  would  be  obtained 
from  the  formula  •  - 


I  AMV  2* 


Here  TJ,  ■  and  u  are-  known  if  the  integral  equations  have  been  solved, 

jL  Q  bV 


and  ^(ru,^)  be  obtained  from  the  numerical  differentiation  formula 


(54)."  Thus  the  pressure  coefficient  is  given -as  a  function  of  x  and  Oi 


-£■— — -  t-  (  V,?  tet/n  -  wf  arc  6 
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+  [(  c&axy  CtnL*6  +  U-tf  AvrSd ]  Aaav*<*  j-  (90) 
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Vrnen  the  motion  is  unsteady,  the  pressure  depends  upon  the  tine 
rate  of  change  of  the  velocity  potential,  so  that  the  instantaneous  state 
.of. ‘motion  does  not  suffice  to  make  the  pressure  determinate.  A- practical 
application,  in  which  the  'Various  potential  flow  solutions  have  been  com¬ 
bined  for  a  body  of  revolution  maneuvering  in  an  arbitrary  trajectory,  is  ■" 
given  in  TM3  Report  No.  987.  [10].  ’  .  • 

When  a  computer  is  not  available,  it  has  been  recommended  that 
•  several  short  .cuts  be  applied,  although  at  some  sacrifice  in  accuracy..'  ’  • 

■  -These  are  illustrated  on  the  same1  body  in  Tables  4,  5  -and  6  for  the  case 
of  axi symmetric  flow  said  compared  with  the  machine- compute d  refeulte  in 
Figure  4.  Also  shown,  in  this  figure  are  the  results  obtained  by  Kaplan's 
-  method  [111*  .  •  • 

•  ,  .  A  Girailor  comparison  is  shown  in -Figure  5- for  the  two-dimensional 

,  form  in  longitudinal  flow.  The  pressure  distribution  obtained  from  the 

values  in  Table  3  are  here  compared  with- the  results  from  Theodorsen's 
method  [12]*  .  ... 

In  order  'to  try  the  procedures  for  a  two-dimensional  profile 
with  a  sharp  trailing  edge,  for  which  the  Kutta  condition  must  be  applied,  . 
a  third  set  of  calculations  were  performed  for  the  Joukowsky  foil  shown 
in  Figure  3*  This  was  derived  from  the  circle  of  radius  a  in  the  $ -plane 
;■  £  •=  a  -  i  +  a.e  10  ' 

•  by  the  transformation  ,  . • 

5  =  [(za-i)C  !?  +  jr  )  -  2<a-i)zJ 

which  gives  a  symmetrical  profile  with  a  sharp  trailing  -edge  at  x  =  -1  and 
a-  rounded  leading  edge  at  x  =  +1.  The  equations  of  the  profile  in  para- 
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The  case  shown  in  Figure  3  is  for  a  =  1.20. 

The  exact  potential -flow  solutions  for  the  foil  are 
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Qx=  (t+ttrt,e)[' o*-20+2  -  a(2-a)  . 

and  the  Kutta  condition  gives 

•  XL  -  _*£!!  (.za.-0  , 

.  •  ’V  a- 

Hence  the  resultant  velocity  \m  the-  body  'is 

and  the  pressure  distribution  is  given  by 

—  _ L.  I  -  (  -£-)*.[  xXn(O-x)  - 

tp(ut-fV'a;  ^  oft/  l  J 
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where  «t  a  arctan  V/U.  •  * 

The  velocities  on  the  Joukowsky  profile  in  longitudinal  flow, 
obtained  by  solving  the  integral  equation  with  the  aid  of  the  IBM  650, 
are  given. in  Table  7.  The  subsequent  calculations  in  this  table  illus¬ 
trate  the  procedure  for  deriving  the  velocity  and  pressure  distribution 
when  the  body  is  .at  an  angle  of  attack;  10  degrees  was  assumed  in  the 
present  case.  The  pressure  distributions  are  compared  in  Figure  6. 
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,  •  Table  2  -  Prelisxinar?  Calculations  .for  IS!  &j0  1’ro'cran 
.'  •' ■  for  r2  =  0.0*i(l~ x*) 


i 

•  ri  .. 

G,  =  — (r2). 
i  2cscx  ' 1 

fiec2y^ 

i 

-0.98940093 

.  .  0.001669079 

•  0.077483091 

4.59697138 

2  . 

-0.94457502 

0.608157598  i 

.  *  i 

0.067421647 

1.55723248 

3 

-0.86563120 

. 

0.017540959 

■  1 

•  0.051390600 

1.15550571 

4 

-0.75540.141 

0 ..026974990  ■ 

0.034484865 

1.04408550 

5 

-0.61787624 

0.034170034 

0.0188709S1 

■  1.01042182 

6  1 

•  -0.45801678 

0.038239704  | 

0.0076855976  . 

1.00154509 

7  . 

•  -0.28160355 

0.039748457 

0.0017865055- 

1.00008029 

8 

•  -0.09501251 

0.039996740 

0.0000686171 

1.00000012 

9  ' 

0.09501251 

0.039996740 

"  -0.0000636171 

1.00000012 

•10 

0.28160355 

0.039748457 

-0.0017865055 

1.00008029 

11 

0.45801678 

*  0.038239704 

-0.0076865976 

1.00154509 

12 

0-61787624 

0.034170034 

-0.018870981 

« 

1.01042182 

13 

.0.75540441 

0.026974990 

-0.-034404865 

1.04408550 

■ 

14 

0.86565120 

0.017540939 

-0.05189,0600- 

• 

1.15350571 

15 

0.94457502 

0.008157598 

■ 

,-0.067421647 

1.55723248 

16 

.  0.98940093 

'  0.001669079 

* 

•  • -Q.077483091 

4.59697138 

of  revolution 


:  Two- dimen si onal 


Body 


■ 

-U..  sec  Y 

Is  * 

V29 

■SB 

Ulssec  y 

xy  • 

B 

•  1.170939 

1.776103 

-1.611283 

1;  3894 27.  a 

-2.659327 

.2 

1.155112 

1.787106. 

-1.551686 

.  1.368001 

-2.72i678 

3 

1.141836 

'  1.801656  ' 

-1.439410 

1. 341782 ' 

-2.783120 

4 

1.118176 

1.825418 

-1.279773  . 

1.301883 

-2.903859 

5 

1.091943 

1.853961 

-1.069772 

1.256147* 

-3.042681 

6 

‘  1.064840 

1.885408 

-0.811441  ’ 

1.209526  • 

.  -3.199096" 

7 

1.041513  ' 

1.914041 

-0.509042  ■ 

1.169843 

-3.340109  . 

8 

1.027651 

1.951756 

-0.173839 

1.146485 

-3.426730 

9 

1.027651 

1.931757 

0.173839 

1.146485 

-3.426730 

10 

1.041513 

1.914041  . 

0.509042 

1.169843 

-3.340110 

11 

1.064840 

1.885409 

0.811441  . 

11209526  * 

,  -3.199096  • 

12  • 

1.091943 

1.853961 

1.069771 

1.256147 

-3.042681 

13  ' 

£.118176 

1.825417 

1.279773 

1.301882 

-2.903859 

14 

1.141836- 

1.801657 

-  1.439411 

1.341782 

-2.783121 

15 

1.155113 

1.787105 

1.551625 

1.368001 

-2.721577 

16 

1.170939 

1.776102 

I.6U283 

' 

1.389427 

-2.659328 

Table  S'  -  Short-cut  procedure  for  axial-flov  pressure  distribution 
on  body  r^  =  0.04(l-x4) 

Calculation  of  matrix  Ik  ^  for  solution  of  linear  equations  by  iteration 
h,j  ='  Rjhtj  +  &:j  ( r't  -  Z'  Hi  kij)  - 


1T3/  V-  ,  . 


1  L  -f 

r.  ; 

fc.. 

OX 

.055669  j 

17.222847 

.125560 

2.669813 

.186290 

.2030571 

•  .233194  ! 

.0340121 

.262805 

.0093735 

.272925 

.0035824 

. 262805  ' 

.0017301 

• .253194 

.0010020 

.186290 

.0006750 

.125580 

.0005187 

.055669 

.0004496 

i  -  . 

R.k.. 

1  Jx  _ 

•958779 
.535275 
.057973 
.007931 
.002463 • 
.000978 
.000455 
.000234 
’.000126 
.000065 
..000025 


i  !  v 


k. . 

R.k. . 

4i  ■ 

•  3  Ji 

4.216571 

.234732 

8.102328 

1.017490 

1.914577 

.556667 

.2605540 

.060760 

.0609713 

-016024 

.0212043 

.  .005787 

.0096816 

.002544 

.0054159 

‘.001263 

-0035702 

.  .000665 

, .0027082 

.000340 

.0025318 

.000150. 

1.017490 


ft? 

$' 
m * 
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Table  5  -  -  continued  . 


r, 

-4 

■  • 

__i_=  5  ... 

. 

1  =  6 

* 

- 

3 

V  - 

VCji 

RAi  ] 

KJi 

ilk.. 

.  3  o* 

R.k.. 

.  111  J* 

- 

■: 

‘j*! 

•  • 

l 

.0996740 

.005549 

•  -.026195 

.04)1048 

.002237 

'  .010907 

i- 

•  2 

.1455404 

.018277 

.038249 

.0531930 

.006681 

-014519 

4  * 

5 

.3144701 

.058583 

'  .002644 

.0921912 

.017174 

.025181 

r 

< 

1 

4 

1.219429 

.  284364 

.320472 

.2523548 

•  .,054184. 

.065415 

j 

5 

5.015236 

1* 517503 

1.317503 

1.058579  ' 

.278200. 

’ . 283913 

> 

\ 

6 

1.055965 

.238199 

.-277513 

5.000000 

1.364625 

1.364625 

Fi¬ 
fe*  _ 

: 

♦ 

7 

.2096917 

.055108 

.055108 

1.058579  • 

.278200 

.288313 

t 

1  .  *  -  * 

8 

-.0759588 

.017242 

.019451 

.2323548 

.054184  • 

.063415 

i- 

.  - 

9 

.0575700  • 

.  .006999 

.009874 

.  .0921912 

-  .017174 

.025161 

r. 

* 

10 

.0246157  ■ 

.003091 

.006469 

.0551980 

. 00G681 

.014519 

1  . 

11 

.0197247 

.001098  ■ 

.005184 

.0401848 

.002237 

.010967 

t 

< 

•  *  t 

i  =  7  . 

• 

» 

i  =  8 

. 

h. 

* 

f 

•r 

; 

' 

J 

/•  V 

.  RJkdi 

.R.k,. 

•  ■  kji  . 

n.k' 

.  3 

Rikji 

K 

K. 

1 

.0197247 

.001098 

.005184 

.0107016 

.000600 

.002514 

f 

2 

.0246157 

.003091 

.006469 

.0129751 

.001629 

.005026 

c  . 

'  3 

.0575700 

.006999 

.009874 

.0183685 

.003422 

.004283 

i 

4 

.0739388 

.017242 

.019431 

'  .0515099 

.007343 

.007348 

l 

} 

% , 

5 

.2096917 

.055108. 

.055108 

.0693567 

.018227 

.016174. 

V 

\ 

% 

< 

6 

1.055965  , 

.280199 

.277513-  j 

.2185553 

.059649 

-0SC965 

? 

i- 

7 

5.013236 

1.317503 

1.317503 

1.185395 

.311002 

.275961 

f 

1  . 

4 

.8 

1.219429 

. 284364 

1520472  !|  5.192052 

1.210755 

1.210755 

t 

V  • 

9 

,  .3144701 

.058583 

■  .082644  !!  1.590250 

.296248 

.370837 

\ 

10 

.1455404 

.018277 

.038249 

.5177928 

.065024 

.120746 

■. 

l 

11 

.0996740 

,.  005549 

,  '  .026195 

.  .3005462 

.016731 

...07.0086 

‘ 

1  ■■ 

.  . 

/ 

1  =  9 


i  =  10 


i 


1 

2 

.3 

4 

5 

6 

7 

8 
9 

10 

11 


k 


ji 


.0056588 

.0066595 

.0090120 

.0142906 

.0274758 

.0680013 

.2421879 

1.446400 

5.909933 

2.329880 

1.057357 


Yji 


.00031.5 

.000836 

.001679 

.003332 

.007221 

.018559 

.063648. 

.337292 

1.1Q0961 

.292586 

.058862 


Vji 


.001054 

.001241 

.001679 

.002662 

.005118 

.012668 

.045117 

.269450 

1.100961 

-.434033 

.196975 


V 


■  .0023318 
.0027082 
.0035702 
.0054159 
.0096816 
.0212043 
.0609713 
.2605548 
1.914577 
8.120328 
4.216571 


Yji 


.000130 
-.000340' 
.000665 
.001263 
.002544 
.005787 
.016024 
.060760 
•  .356667 
i  1.017490 
J  .234732 


Rikdi 


.000293 
.-000340 
.000448 
.000680 
.001216 
.002663 
;007657 
.052720 
. 240433 
1.017490  j 
.529517  J 


•K 


K 
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l’abLc  G  -  Short-cut  yrvcjfuiro  *Vr  r  *=  0.04  (I*-;:  ) 
Solution  Ojl  equations;  pressure  uisoribution 


3 1  -  4  co°  ^ 


cn+l,i  3  Gni 


+  Jj  i  I 

m7 


f  <2  •  ta  ‘u 


E1J 

C2j 

r,» 

C-, 

oj 

E3j 

*Di 

JtL  W 

.01286 

1.15904 

.00164 

1.16068 

-.00026 

1.16042  0.49445 

.01312 

1.14864 

.00304 

1.15168 

.00092 

1.15260  -.10275 

.00690 

1.11500 

.00118 

1.11618 

.00010 

1.11628  -.20526, 

.00084 

1.07690 

.00014 

1.07704 

.00004 

1.07708  -.15620 

-.00462 

1.04262 

-.00088 

1.04174 

-.00016 

1.04158  -.08482 

-.00692 

1.02816 

-.00132 

1.02684 

-.00024 

1.02660  -.05391 

-.00462 

1.04262 

-.00088 

1.04174 

-.00016 

1.04158  -.08482 

.00084 

1.07690 

.00014 

1.07704 

.00004 

1.07708  -.15620 

.00690 

1.11500 

.00118 

1.11618 

.GC010 

1.11628  -.20526 

.01312 

1.14864 

.00304 

1.15168 

.00092 

1.15260  -.10275 

.01286 

1.15904 

.00164 

1.16060 

-.00029 

1.16042  .49445 

J  Gl<  **  I 
u  J 

1  1.146180 

2  1.135520 

3  1.108104 

4  1.076066 

5  1.047236 

6  1.035072 

7  1.047236 

8  1.076066 

9  1.108104 

10  1.135520 

11  1.146100 


Tabic  7  -  Velocity  and  pressure  distributions  for  a  Jouicowoky  foil 
!  I  ..  //  f  t.2\  I  r-T7T7T-| 


computed 


see  f 
'"exact 


.9894 

•000o 

i8495 

.8565> 

.9446 

.0031 

.8504 

.0519, 

.8656 

.0114 

.0713 

.8706 

.7554 

.0272 

.9108 

.9163 

.6179 

.0511 

.9574 

.9639 

.4580 

.0820 

1.0163 

1.0194 

.2016 

.1172 

1.0015 

1.0805 

.0950 

.1526 

1.1467 

1.1441- 

.0950 

.1836 

1.2008 

1.2071* 

.2016 

.2056 

1.2670 

1.2666 

.4580 

.2149 

1.3200 

1.3202 

.6179 

.2090 

1.3662 

1.3662' 

.7554 

.1874 

1.4059 

1.4038 

.8656 

•  .1509 

1.4313 

1.4326 

.9446 

.1026 

1.4503 

1.4523 

.9894 

.0463 

. 

1.4651 

1.4652' 

1 

p KW) 

ft  *  o 

cd,$utcd[ '“criict" 

.2810  .3006 

.2632  .2792 

.2404  .2599 

.1911  .1812 

.1138  .1018 

.0064  '  .0002 
-.1258  -.1257 
-.2738  ;  -  .2681 
-.4316  -.4277 
-.5929  -.5918 
-.7422  -.7426 
-.8508  -.8508 
-.8675  -.8673 
-.7011  -.7042 
-.2160  -.2194 
.6189  .6199 


n pit  p/(W> 

iOn) /WiFT0,nputc1d><tM  10 

»  -l  r  '“i^pper  lower 


.0090 

.0472 

.1157 

.2145 

.5427 

.5004 

.6855 

.8934 

1.1174 

1.3483 

1.5765 

1.7913 

1.9818 

2.1381 

2.2518 

2.3171 


.0621 

.1438 

.2285 

.3181 

.4150 

.5226 

.6451 

.7879 

.9589 

1.1704 

1,4440 

1.8189 

2.3765 

3.3146 

5.2762 

12.1974 


.3179 

.3257 

.3213 

.3010 

.2617 

.2058 

.1424 

.0840 

.0413 

.0269 

.0612 

.1717 

.3887 

.7151 

.9943 

.5106 


.2873 

.2538 

.2027 

.1250 

.0102 

-.1494 

-.3544 

-.6025 

-.8977 

-1.2483 

-1.6596 

-2.1310 

-2.6471 

-3.1418 

-3.5946 

-2.6690 
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■  APPENDIX  1  ■ 

Formulas  for  Numerical  .Differentiation  and  Interpolation 


let  us  suppose  that  a  Gauss  quadrature  formula  of  order  II 
las  "been  used  to  solve  an  integral  equation.  The  solution,  will  consist 
of  a  set  of  numbers  g(x^)'  where  the  x  axe  the  zeros  of  the  Legendre 
polynomial  of  degree  H.  It  ii:  necessary,  occasionally,  to  differenti¬ 
ate  the  product  of  g(x)  by  Some  given  function.  Since'  it  is -desired 
to  maintain  a  completely  arithmetical  procedure,  a  formula  for  the  - 
values  of  the  derivative  of  this  pr6<luct  ,at  the  point  x^  will  be  de¬ 
rived.  ■  .  ■  ■  ,  -  . 


■Denote  the  product  function  b>  h(x)  and  put  h^  -  h( x^ } .  The 


U  numbers'll^,  h0  —  h;.  are  presumed  to  be  known.  Lagrange's  iaterpo-1 


latioa  formula  [13]  for  the  polynomial  through  the  points  (x.,  b  )  is 

.  0  U 


fu.. 


■  h  Cx)  =  PJ*)  Z  - - 

"  J=' 


t  ty) 


(91) 


where  the  prime  denotes  differentiation -with  respect  to  the  argument. 
*  . 

Differentiating  (91)  yields  ' 

N 


■  u’m  -  z  h, 

or,  putting  x  =  x^  and  h^  =  h,(x.. ),  we  obtain 


*>'  =  /S,  a*ihi 


(92) 


where 


a.,  - 


lJ  .  (*4-Xj)pK^> 


ift  j 


(93) 


S3 


m 


The  interpolation  formula  (91)  also  furnishes  an  ' 
arithmetical  praee&m:  Snie-  ewulmtin®;  the'  •Integrals  in  (ff). 
We  obtain . 


The  coefficients  C. .  have  the  properties 


'lj  +  l+l,  #-^+l~  x  ~ 2~ 


vhere  the  Rj  are  the  weighting  factors  in  Causs’s  quadrature  formula 
(SB)/  Also  we  h.ave  from  (99)  •  .  ‘ 

f '  h(X)  j*  z=.  J  h(*)J*  -  J  h  lf)  dx 

**  ■  -  •  -»  •;  .  •  •  (100). 

N  AT 

=  (  ni  ~  C‘i  ~  2-t 


Values  of  the  matrix  for  N  =  16  are  given  in  Table  9  for  1  =  1,  2, 
•••  8,  and  j  ■  1,  2,  •••  IGj  The  values  for  i  =  -9,  10,'  •••  16.may  be 
obtained  from  the  tabulated  values  by  means  of  (99).  An  alternative  , 
procedure,  which  avoids  the  necessity  of  'computing  these  additional 
numbers,  is  the  following  one: 

f'  h(r)<h  =  H-iti)  hi  *  (101) 

It  is  seen  that  Et^s*  (100)  furnish  similar  expressions  for  evaluating 
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Anpcndix  2 

.  Program  for  the  IBM  650  Automatic  Computer 

by 

Matilde  thcagno 

The  IBM  650  io  a  computer  of  moderate  speed  and  a  storage 
capacity  of  2000  words  which  may  Gerve  as  numbers  or  inatructiono. 
Because  of  its  limited  capacity,  each  of  the  following  two  programs 
was  subdivided  into  two  parts.  The  continuity  of  a  program  is  not 
disturbed,  however,  since  the  second  part  is  read  in  automatically 
after  the  first  is  finished.  The  programs  are  written  in  the  form 
of  FORTRAN  statements,  a  semi-mathematical  language,  which  are  trans¬ 
lated  by  the  computer  into  a  imchine-language  program. 

The  first  program  solves  the  three  integral  equations  for 
axial,  transverse  and  rotational  motions  of  a  body  of  revolution. 

The  required  input  data  are  the  dopes  and  ordinates  of  the  body  at 
16  specified  points.  In  the  first  port  of  this  program  a  "corrected" 
matrix  for  axial  flaw  is  comimtcd  and  various  combinations  of  the 
elements  of  this  nutri;:,  needed  for  all  three  problems,  axe  obtained. 
In  the  second  part  a  system  of  linear  equations,  having  the  matrix 
computed  in  the  first  part,  is  solved  by  iteration  to  give  the  solu¬ 
tion  for  axial  flow.  Then  a  second  matrix,  which  serves  for  both  the 
transverse  and  rotational  cases,  is  derived  and  used  to  obtain  solu¬ 
tions  (also  by  solving  linear  equations  by  iteration)  for  the  latter 


two  cases 


The  second  program  solver,  the  two  integral  equations  i'o>- 
longitudinal  and  rotational  motion  of  a  symmetric  two- dimensional 
form.  Again  the  first  part  of  the  program  derives  only  quantities 
associated  with  the  matrix  for  longitudinal  motion.  In  the  second 
part  the  matrix  for  rotational  motion  is  obtained  and  solutions  for 
the  two  cases  are  computed  by  solving  seto  of  linear  equations  by 
iteration. 

A  "Pause"  statement  is  included  to  pormit  the  program  to 
be  stopped  after  a  solution  for  a  particular  mode  of  notion  has  been 
obtained.  Two  subroutines  are  required,  the  square-root  subroutine 
for  both  programs  and  the  natural- log  subroutine  for  the  body  of 
revolution. 

The  times  required  for  processing  the  programs  on  the  IBM 
650  at  the  State  Univorsity  of  Iowa  (which  is  presently  not  equipped 
with  "floating  point")  are  as  follows: 

Body  of  revolution  (3  integral  equations)  100  minutes 

Two-dimensional  form  (2  integral  equations)  75  minutes 


VlMtHSlW  C( 16,1s),  S(16),  Z(16),  X(l6),  Y(l6),  R(l6),  W(16), 
H(l6),  Si(l6),  S2(l6),  A(16),  B(l6),  T(l6) 

1  HEAD,  X,  Y,  R 
DO  3  I  <=  1,16 

2  S(I)  =  1.  -  X(X)  **2 

3-  H(I)  .  Y(l)/s(l)  ,  _ 

.a  !  ■  DO  6  I  =  1,16  ' 

-•  l  . 

4  50  6  J  =  1,16  - 

5  HE  =  S(J)  *  11(1) 

D  »  (X(J)  -  X(I))  **  2 

6  C(I,J)  =  R(j)*  YE/SQKEF((D  +  YE)  *r3) 

DO  8  I. «  1,1$  . 

7  B  a  SQRH'(l.  -  H( I ) ) 

•U  a.LKF((l.  +B)/(1.  -  E)) 


A(l)  =  ((2.*E  -  D*Il(l))/Ex*3 

8  B(l)  =  H(I)*(6.*E  -  (2.  +  1I(I))*U)/B**S 
DO  11  I  a  1,16' 

9  Sl(l)  =  0.0 
S2(I)  =  0.0 

DO  11.  J  a  1,16 

10  V  a  C(l,j)  • 

S1(I)  a  S1(I)  +  X(l)*B(l)/l6.  -  V*(X(j)  -  X(l)) 

11  S2(I)  a  S2(l)  +  A(I)/16.  -  V 
DO  14  I  a  1,16 


-J.  JV  i  jVAT^^..l^UJ%^>»xr^U^U^ZyAt5^lWiJy^rK.\.  X.UX.I^ulI>UiriL!3  A1.-.X.  w  jO^a  l^HXrJj.T1  Jy1 


kJ 


12  DO  14  J  =  1,16  .  -  ■  \ 

iJ5>  ED  «:  (X(j)  -  X(3)>»2  +T(J) 

.14  C(I,J)  =  H{j)*l(j)/SQRTF(i>**3) 
DO  18.  I  ■  1,16 
/  15  S(l)  *  0.0 
Z(l)  =>  0.0 

.  .  DO  17  J  -  1,16  ‘  •  . 

16  .  V  =  C(I,J)  . 

.  Z(I)  «  2(1)  +  V*(X(j)  -,x(l)) 

17  S(l)  =  S(I)  +7  .’ 

S(l)  -  S(l)‘  +  S2(I) 

H('l)  =  Z(I)  +  S1(I) 

18  C{l,l)  =  C(I,I)‘  +  S2(I)  ’  .  • 

..  D02-  I  *•■  1,16  -  .' 

19  READ,  W'  -  .  .  .'  • 

DO  20  J  =  1,16  ' 

20  C(I,J)  =  C(I,.J)  +  W(j)*Sl(l)  • 
21'  Z(I)  =.  2.- 


Body"  of  Revolution  -  Part  2 


■DIMENSION  -  Same  as  in  Part  1 


■  M  =  1 


24  DO  26  I  *  1,16 

25  A(I)  =  Z(I)/S(I) 
AP  =  A(l)  ‘ 


-irvt*  -ji*  -2- 

•  y  v 7 -v^y n*r>—  >—  ,•  » - 


r  ^  ^  .•  /  ^  - '  * 


JU3£&KiCIUUL2U  AKAUiUiiJLUJUi'  TLKK^^WMXVXkTvlZX  kTCflfTCESVVi 


i 

i 

f,  * 


•  ( 


26  fUHCH,  I,  AP 

*  k  =  i 

27  'DO  51  I  -1,16 

28  S1(I)  «*  0.0  ’ 

DO  50  J~=.  1,16  -‘ 

29.  D  =  Z(l)/l6  =■  ’" 

30  S1(I).=  Sl(l)  +  D  -  C(I>J)*A(J)  .  • 

E  =  Sl(l)/S(l). 

B(l)  »  A(I)  +  E  ■  : 

bp  -  B(i)-  ,  ;  ,  .  ■  ' 

31  PUNCH,  I,  BP,  E  ' 

DO  55  1  =  1,16  •  - 

32  S1(I)  *  0.0  '  1  . 

.  DO  54-  J  =  1,16  ’  .  • 

33  D:  Z(l)/l6 

•  34  S1(I)  =  Sl( -I )  +  D  -  C(I,J)*£(J}  ' 

E  *  si(i)/'s(i)  . 

A(l)  =  B(I)  +  E‘  • 

AP  =  A(I)  ' 

35  PUNCH,  1/  AP,E. 

K  a  K  +  1 , 

.  IF(K  -  6)  27,  36,  56 

36  CONTINUE  - 
.  M  = M  +1 

PAUSE  - 

I?(K  r  3)  50,  51,  53  ■ 


.50  BEAD,  A 

DO  39  I  =  1>16 


.57  '  S2(I)  «  O'.O 

DO  30  J  a  1,16  , 

38  S2(I)  =  S2(I)  +  C(l,j)*A(j) 

■  K(I)  =  H(I)  -  62(1)  *  ■ 

S(I)  «  S(I)  -  4 

39 ;  T(I)  «  4.*X(l)  +  2.*S2(I)  ■  " 

•  •  -  DO  42  I  . =  1,16  ■  .  ’ .  ;  .  ■  • 

40  DO- 42  J  i  1,16  ;  ’ 

;  .41.  D  =  (X(J)  -  X(I))^2  +  Y(j)  ■ 

E  =  D  -  3.*  Y(J)  +  3.  (X(J)  ~  X(.l))*A(j)  '  . 
G  =  R(J>  Y(j)/SfiftTF(D#*5)  •  • 

42  C(J,J)  =  BkG 

DO  46  1  =  1,16  ■* 

•  43  Sl(l)..=  0.6  .  '  .  ■ 

■S2(I)  =  0.0'.'. 

•  DO  45  -J.  =  1,16  '  •’ 

-  44  V  =  C(I,J) 

S1(.I)  =  S1(I)  +'K(l)/l6.  -  V*(X(J)  -  X(I)) 

45  S2(I)  =  S2(I)  +  S(l)/16.  -  V 

46  C(l,l)  =  0(1,1)  +  S2(l)  *;  ■  ’  ' 

DO  49  I  =.1,16 

47  READ,  W  '  •  ; 

DO  48  J  =  1,16 

48  C(I,J)  =  C(I,J)  +  V(j)*si(l) 

49  Z(I).  =  4...  • 

GO  TO  24  . :  -  •  .  '  ' 

51  '  DO  52  I  »  1,16 


i- 


-sa-  ■ 

52  2(l)  ««  T(l)  ’  .  •-  ' 

00  TO  24  ‘  .  ‘ 

53  CONTINUE  •  -  -  .  ■•  •  • 

STOP  ‘  ’  ■  .- 

1  •  Potential  Floy  Afcout  A  Symmetric  Tvo^Dlmcnsional  Form  ,  . 

Part  1  .  -  '  \  ' 

DIMENSION  "  C(l6,16) ,  S(l6),  Z(l6),  2(16),  2(16),  R(l6),. W(16),  H(l6), 
Sl<16),  S2(16),  A(16),  B(16)  ..  /\  .  ‘  *" 

•  1  -READ,  X,  Y,-  R  ' '  *  .  v  ' 

P  *  3.1415927  •  "  - 

DO  3  I  =  1,16  "  .  .  ■ 

2  "  S(r)  =>  1.  -  X(l)**2  ■  ‘  • 

3  H(I)  *  (Y(l)**2)/s(l)  ....  .  . 

DO  6  I  *  1,16  ■  •  ' 

4  DO  6  J  *  1,16  .  ‘  '  *  ** 

■5  ■  YE'  =  S(j)*H(l)‘.  ‘  ■  * 

D  *  (X(j)  -  X(l))#*2  ‘ 

.  ‘  ,  '  «  \ 

6  C(I,J)  »  R(J)*S«RTF(YE)/(D  +  YE) 

DO  8  I  *  1,16  • 

7  U  =  SQRTF(H(I)) 

A(l)  a  P/(l.  +  U)  •'  *  .  '  ’  ‘ 

8  B(l)  -  -'P*U/(l.  +  U)**2  . 

DO  11  1  *1,16 

9  S1(I).=  0.0  '  ...  . 

S2(I)  =  0.0.  .  . 

D0  11  J*  1,16 


..  jo  t«  c(r,j)  ■  ' 

Sl£c)  *  Sl(E]i  +•  X(lUB(l)/lS.  r  V*(fcCJj-  -  *(,£}$ 
n  S2(i>  *  ^(i)  +  a(i)/ig.  -  v 
'  DO  14  I  rn  1,16 

12  DO  14  J  -  1,16  '  '  •  ’ 

•  •  ’13  D  .*  ’  (X(J)  ~  X(l))#*2  +  Y(J)**Z 

,14  C(I,J)  »  R(j)*Y(J)/D  ■ 

DO  18  I  -  1,16  •  ,  \ 


.15  Z(l)  «  0.0. 

S(I)  -  0.0  ■  • 

•  i  .  •  .  •'u 

•  'DO  18  J  -  1,16  .  . 


i  16  V  ».C(I,j) 

* 

,-:  \  .  •  .  . .  z(i)  *  z(i) +v*(x(J)  -  x(i))  .: 


■17  S(I)  *  S(I)  +  7 

S(I)  ■«  S(I)  +>S2(l) 

18  H(l)  =  Z(l)  +  S1(I) 

’  DO  21  I  s  1,16  ’  : 

19  -READ,  W  •■  ■ 

DO  20  J  »  1,16 

20  C(I,J)  »  C(I,J)  +  W(J)*S1(I)- 

-  C(.I,I)  »  C(1,I)  +  S2(I) 

■  21  z(i)  a  p  • 

STOP  •  •’ 


'  TW-Pin^tiowai.  form  -  fttrt  ft 

XWtfft&tffi  (C^tSS,!®),  S(M3),  '<0[DiS),  E(ES)s  B(1S>),  WftnffiJ 

11(16),  S1(16),  S2(l&),  A(16),  B(16) 

24  BO  26  I  =  1,16  *  • 

25  A(l)  =  Z(X)/s(l)  •  ..  .:  ■ 

AP  =  A(I)  ‘  I  ‘  ■ 

26  PUJ1CH,  I,  AP 

•  K  .=  1  ,  ■'  ‘  ’  ...  . 

27  BO  51  I  =  1,16  ■  • 

28  S1(I)'=0.0  .  .  ’  '  . 

DO  30  J  =  1,16 

29  S*  Z(l)/l6  .  ' 

-30  S1(I)  =  S1(I)  .+  D  -  P{I,J)*A(J)  . 

,3»$1(I)/S(I)  • 

.  .  B(.I)  =  A(l)  +  E  *  *  . 

BP  =  B(I)  '  . 

31  PUHCH,  I,  BP,  E  ■ 

.  DO  35  X  =  1,16 

32  S1(I)  =  0.0 
■  DO  34  J  =  1,16 

•  w  T  *  • 

33  D  =  Z(l)/l€.  . 

34  S1(I).  «  S1(I)  +  D  -  C(I,J)*B(J)  . 

E  =.  S1(I)/S(I)  ' 

A(I)  =  B(I)  +  E  •  ' 

• .  AP  =  A(I) 

35  PUHCH,  I,  AP,  S 
K  =  K  +  1 

ZF(K  -  6)27,  36,  36  - 
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